Orbital Ordering and Spin-Ladder Formation in La 2 RuOs 
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The semiconductor-semiconductor transition of La2RuOs is studied by means of augmented spher- 
ical wave (ASW) electronic structure calculations as based on density functional theory and the local 
density approximation. This transition has lately been reported to lead to orbital ordering and a 
quenching of the local spin magnetic moment. Our results hint towards an orbital ordering scenario 
which, markedly different from the previously proposed scheme, preserves the local S — 1 moment 
at the Ru sites in the low-temperature phase. The unusual magnetic behaviour is interpreted by 
the formation of spin-ladders, which result from the structural changes occurring at the transition 
and are characterized by antiferromagnetic coupling along the rungs. 

PACS numbers: 71.20.-b, 71.30.+h, 71.70.Ch, 71.70.Gm 

Keywords: density functional theory, orbital ordering, spin-Peierls transition 



The orbital degeneracy of d-shell atoms lays ground 
for numerous exciting phenomena observed in transition- 
metal compounds Orbital fluctuations as well as or- 
bital ordering lead to extraordinary ground states, low- 
energy excitations and phase transitions. Well known ex- 
amples for such ordering phenomena are the perovskite- 
based manganites 0, Q , and the antiferro-orbital struc- 
ture in KCuF 3 Q. Increased complexity is observed 
for compounds, where the orbitals couple to the spin or 
charge degrees of freedom. This situation has been stud- 
ied by Kugel and Khomskii for magnetic systems [f| . Or- 
bital and magnetic ordering has also been found in the 
triangular chain magnet Ca3Co20g, which is character- 
ized by an alternation of low-spin and high-spin sites [(| . 
In contrast, interplay of charge and orbital order as well 
as singlet formation has been demonstrated to play a sig- 
nificant role in the Magneli phase T14O7 0. While or- 
bital ordering has been primariliy studied for transition- 
metal compounds of the 3d series, interest in the 4d ox- 
ides has considerably grown. Prominent examples for 
such oxides are Ca2Ru04 8], Sr 2 Ru04, which shows 
superconductivity below T c » 1.5 KW, and SrRu03, 
which is ferromagnetic below 160 K [lOj. 

Recently, focus has centered on the new ruthenate 
L^RuOk, which shows a first-order phase transition near 
160 K |ll|. This semiconductor-semiconductor transition 
is associated with a slight increase of the band gap from 
« 0.15 eV to about 0.21 eV. In addition, it is accom- 
panied by strong changes in the magnetic properties as 
well as by a transformation from a monoclinic to a tri- 
clinic lattice. High-temperature Curie- Weiss behaviour, 
with /i c ff = 2.53/zb and O = — 71 K, is attributed to 
the low-spin (5* = 1) moments of the Ru 4+ ions. At 
the transition, the susceptibility drops to a small, nearly 
temperature-independent value with a slight upturn at 
lowest temperatures assigned to free intrinsic spins or ex- 
trinsic impurities 0,^2- Furthermore, from the absence 
of any field dependence of the magnetic susceptibility for 
fields up to 9 T a complete quenching of the magnetic mo- 
ments was deduced and attributed to an orbital ordering 
of the Ru ions. 




FIG. 1: Crystal structure of ht-La2RuOs viewed along the in- 
direction. La, Ru, and O atoms are displayed as big, medium, 
and small circles, respectively. 



A schematic representation of the crystal structure is 
displayed in Fig. ^ Arrows indicate the random orien- 
tation of the spin moments in the high-temperature (ht- 
) phase. The structure is built from octahedral chains 
parallel to the c-axis, which form a zigzag-like pattern 
in the afr-plane, thus giving rise to double-layer slabs, 
which are separated by LaO-layers [Til IT3L ll^ . The 
main features of the crystal structure are preserved in the 
low-temperature (It) phase [liT Il4|. However, according 
to the neutron data the magnetic transition is accom- 
panied by pronounced local structural changes leading 
to alternating shortenings and elongations of the Ru-O- 
Ru bond lengths Q. While the Ru-Ru nearest neigh- 
bour distances are quite similar to each other in the ht- 
phase, the lt-phase is characterized by short and long 
distances alternating both within the plane and along 
the c-axis. The in-plane octahedral Ru-0 bonds, which 
range from 1.94 to 2.06 A in the ht-phase, fall into two 
short and two long bonds ranging from 1.87 to 1.97A 
and 1.98 to 2.10A, respectively. The bonds parallel to 
the c-axis are also larger than 2.0 A. In the rotated co- 
ordinate system sketched in Fig. ^ the long and short 
bonds evolving in the lt-phase are along the x- and y- 
axis, respectively. The structural changes lead to two 
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inequivalent Ru sites, which alternate along the c-axis as 
well as the zigzag-like in-plane pattern. ^From structural 
considerations, Khalifah et al. concluded that in the lt- 
phase the Ru Ad yz orbitals are depopulated, leading to 
the configuration d xy d xz d yz , hence, an S — state. In 
contrast, in the ht-phase the t2g states would be nearly 
degenerate and Hund's rule coupling gives rise to the 
S = 1 electron configuration of d xy d xz d yz , d xy d xz d yz , or 
d xy d xz dy Z . Thus, the phase transition was interpreted as 
an orbital ordering transition with a complete loss of the 
local magnetic moment. Finally, the occurrence of an 
inelastic peak at about 40 meV as observed in neutron 
scattering was assigned to the formation of a spin gap in 
the lt-phase Underlining the importance of inter- 

site interactions, Osborn pointed out that the observed 
inelastic response should be rather attributed to singlet 
to triplet type excitations |15| . This point of view was 
supported by Khomskii and Mizokawa 

Concentrating on the above mentioned discrepancies in 
interpreting the data for the lt-phase, we report on den- 
sity functional calculations as based on crystal structure 
data for both phases Q. As expected, we find strong 
changes of the orbital occupations coming with the struc- 
tural transformation. However, in contrast to the pro- 
posal by Khalifah et al., these changes conserve the local 
5 = 1 moment. From spin-polarized calculations it is 
inferred that the suppression of the susceptibility in the 
lt-phase results from a spin-Peierls-like transition coming 
with the formation of spin-ladders with antiferromagnetic 
coupling on the rungs. 

The calculations were performed using the scalar- 
relativistic augmented spherical wave (ASW) method 
fvA ITif . The large voids of the open crystal structure 
were accounted for by additional augmentation spheres, 
which were automatically generated by the sphere geom- 
etry optimization (SGO) algorithm 19]. The Brillouin 
zone sampling was done using an increased number of 
up to 1024 and 2048 k-points in the irreducible wedge of 
the monoclinic and triclinic Brillouin zone, respectively. 
While previous calculations (Ref. |2(J]) were based on 
the local density approximation (LDA) , the present work 
used additionally the generalized gradient approximation 
(GGA) [21| as well as a new version of the ASW code, 
which takes the non-spherical contributions to the charge 
density inside the atomic spheres into account [22|. 

Calculated partial densities of states (DOS) are dis- 
played in Fig. |21 for the high-temperature (ht) and low- 
temperature (It) structure. Contributions from both Ru 
Ad and O 2p states are included. For the oxygen contri- 
butions we have distinguished the atoms 0(1) to 0(4), 
which form the RuOg octahedra, from those of atom 
0(5), which is at the center of the La40 tetrahedra. Not 
shown are the La 5d and 4/ states, which give rise to a 
sharp peak at about 5 eV and smaller contributions be- 
tween —4.5 and — 1.5eV. All other orbitals play only a 
negligible role in the energy interval shown. 

In Fig. |3 the lowest group of bands extending from 
—6.9 to — 1.5eV derives mainly from the O 2p states. In 
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FIG. 2: Partial DOS as resulting from the ht- (top) and lt- 
(bottom) crystal structure. 



contrast, the Ru 4c? states are found mainly in the energy 
interval from —1.1 to 0.5 eV as well as above 1.5 eV. Due 
to the octahedral coordination of Ru by oxygen atoms 
these states are split into t2 g bands near the Fermi energy 
and e g states well above. Hybridization with the O 2p 
states is extremely strong and leads to large p-admixtures 
to the DOS near Ep, which are of the order of « 40% of 
the Ru Ad contributions in that interval. 

The partial DOS of the 0(5) 2p states deviates con- 
siderably from those of the remaining four oxygen atoms 
and shows a steady increase between —4.5 and — 1.5eV 
as well as a sharp drop at the upper edge. La 5c! and 
4/ admixtures in this interval of the order of 15% of 
the 0(5) 2p contribution are attributed to strong hy- 
bridization coming with the formation of the La40(5) 
tetrahedra. The large band width of these rather local- 
ized tetrahedral states reflects their considerable extent 
in space. 

In general, the similarity of the crystal structures of 
the ht- and lt-phase is well reflected by the similarity of 
the partial DOS. Yet, distinct changes are observed for 
the Ru Ad ti g states. In particular, this group of bands 
becomes broader and the strong peak at Ep, seen for the 
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FIG. 3: Integrated partial Ru 4d t2 g DOS as resulting from the 
ht- (top) and It- (bottom) crystal structure. Results for the 
Ru(l) atom of the lt-phase are very similar to those for atom 
Ru(2) as shown here. Orbitals refer to the rotated coordinate 
system depicted in Fig. 

ht-structure, turns into a pronounced dip. Actually, as 
can be seen on closer inspection of the band structure, a 
tiny gap is opened in the lt-phase. 

The differences between the electronic properties of 
the ht- and lt-phase become much clearer from a de- 
tailed analysis of the near-Ep states. For the notation 
of the orbitals the rotated coordinate indicated in Fig. 
is used. In this system, the partial densities of states of 
all three Ru t 2 g bands display the expected rather similar 
behaviour for the ht-phase. As a consequence, integra- 
tion these partial DOS leads to almost identical orbital 
occupations as shown in the upper panel of Fig. Note 
that the t2 g occupations do not vary from to 2 in the 
energy interval shown as would be expected from pure d 
states. Due to the strong p-admixture to these bands the 
energy variation of the orbital occupations appears to be 
somewhat reduced and eventually has to be translated to 
the ideal picture of pure d-states. Obviously, filling the 
nearly degenerate t 2g bands according to Hund's rules 
leads to one of the electronic configurations d xy d xz d yz7 
dlydlzd 1 ^, or d\ y d\ z dl z as proposed in Ref. 



A noticeably different situation is obtained for the It- 
structure. Integrated partial densities of states for the 
atom Ru(2) are shown in the lower panel of Fig. [3J The 
corresponding curves for atom Ru(l) are very similar and 
thus not displayed. Strong similarities between the d xy 
and d yz partial DOS are observed. In contrast, the d xz 
state deviate substantially. In particular, occupation of 
this orbital is much larger as compared to that of the 
other states as is expected from the elongation of the 
Ru-O-bonds along the local a;-axis. Transferring these 
findings to the ideal picture of pure d states leads to the 
electronic configuration d xy d xz d yz . This is in strong con- 
tradiction to the d xy d xz d yz state proposed by Khalifah et 
at, who claim the full occupancy of the d xy orbital at the 
expense of the d yz state. Yet, the similar occupations of 
the d xy and d yz orbitals obtained from the present calcu- 
lations are more in line with bond-length considerations 
using the Ru-0 distances discussed above. Our results 
have important consequences for the magnetic moments. 
While Khalifah et al. propose the complete quenching of 
the local moment, the present calculations clearly reveal 
the conservation of the S = 1-spin moment through the 
phase transition. However, since the orbital degeneracy 
has been lifted by the structural distortion coming with 
the triclinic phase, this moment is carried exclusively by 
the d xy and d yz orbitals. 

The identification of an S = 1-state for the triclinic 
structure motivated additional spin-polarized calcula- 
tions for the lt-phase. Long-range ferromagnetic order 
can be ruled out from the low-temperature susceptibil- 
ity data. In addition, the fact that no extra reflections 
were detected in the neutron diffraction data for the lt- 
phase allows only for antiferromagnetic order with oppo- 
site moments at the inequivalent Ru sites. This situation 
suggests a singlet ground state which, with the discussed 
structural and orbital transitions, may be realized as a 
spin-Pcierls-like state. We have simulated this state by 
starting from opposite moments at the Ru(l) and Ru(2) 
sites. Indeed, our calculations resulted in a self-consistent 
solution, which comes with an energy lowering of 5 meV 
per Ru pair as compared to the spin-degenerate case. In 
addition, an optical band gap of 0.17eV is obtained in 
LDA. GGA calculations result in 0.20 eV, which is very 
close to the experimental value. 

The calculated magnetic moments arise to equal parts 
from the d xy and d yz orbitals thus confirming the expec- 
tation from the previous spin-degenerate calculations. In 
total, local moments of 0.77[1b and —0.73/j.b are obtained 
at the Ru(l) and Ru(2) site, respectively. Together with 
small contributions from the nearest neighbour oxygen 
atoms a magnetic moment of ±0.85/zs per octahedron 
results, which increases to ±1.06/is per octahedron on 
going from LDA to GGA. Interestingly, despite the ab- 
sence of any symmetry constraint, both types of calcu- 
lations lead to a compensation of the magnetic moments 
at neighbouring octahedra and, hence, to an exactly van- 
ishing magnetic moment per unit cell. This finding is in 
excellent agreement with the suppression of the magnetic 
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FIG. 4: Crystal and spin structure of lt-La2Ru0s viewed 
along the c-direction. La, Ru, and atoms are displayed 
as big, medium, and small circles, respectively. 

susceptibility below the phase transition. 

Finally, we performed spin-polarized calculations also 
for the ht-phase, assuming the same type of ordering 
as for the lt-structure. As a result, a metallic solution 
with magnetic moments of ±0.76/is per Ru atom was 
obtained, however, with a total energy much higher than 
that of the corresponding spin-degenerate solution. Tak- 
ing into account both the GGA and the non-spherical 
contributions we arrived at moments of ±1.11/zb per Ru 
atom and an optical band gap of 0.02 eV. The ener- 
getical instability of this solution agrees with the fact 
that no long-range magnetic order is observed for the ht- 



phase. The phase transition to the lt-phase may thus 
be regarded as a spin-Peierls-like transition. Since the 
octahedra form chains parallel to the c-axis, the struc- 
tural transformation leads to the formation of spin lad- 
ders with antiferromagnetic coupling along the rungs as 
indicated in Fig. 01 

To conclude, electronic structure calculations for the 
ht- and lt-structure of L^RuOs reveal strong orbital 
ordering for the latter. While Hund's rule coupling 
within the degenerate ti g manifold leads to a d xy d xz dy Z , 
d\ y d xz dy Z , or d\. y dl. z dt z state in the ht-phase, the struc- 
tural changes associated with the triclinic structure cause 
substantial orbital ordering and drive the system into a 
d xy d xz d yz configuration. While well reflecting the Ru-0 
bond lengths the latter preserves the local S — 1-moment 
and thus is in contradiction with previous interpretations 
of the lt-phase. Spin-polarized calculations for the low- 
temperature phase reveal compensation of the local mo- 
ments due to their antiparallel alignment along the short 
Ru-Ru bonds. As the formation of short and long in- 
plane bonds in the triclinic phase leads to an effective 
pairing of chains, a spin-ladder system is generated. 
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